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We derive the Ginzburg-Landau theory of unconventional singlet superconductors in the presence 
of a Zeeman field and impurities, to examine the resulting Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) 
phases. We show that the behavior of the FFLO phases in unconventional superconductors in the 
presence of impurities is qualitatively different from that found for .s-wave superconductors. 
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I. INTRODUCTION 

In 1964 Fulde and FerrelH. and Larkin and Ovchinnikovi demonstrated that a superconducting state with an order 
parameter that oscillates spatially may be stabilized by a large applied magnetic field or an internal exchange fielcL 
Such a Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) state was subsequently shown to be readily destroyed by impuritiesEl 
and has never been observed in conventional low-T c superconductors. The question of observing an FFLO phase in 
an unconventional superconductor has only been addressed more recent ly n la, particular organic, heavy fermion, and 
high T c superconductors appear to be promising candidates for such statescftJ. These new classes of superconductors 
are believed to provide conditions that are favorable to the formation of FFLO state, because many of them are (i) 
strongly type II superconductors so that the upper critical field H C 2 can easily approach the Pauli paramagnetic limit; 
and (ii) layered compounds so that when a magnetic field is applied parallel to the conducting plane, the orbital 
effect is minimal, and the Zeeman effect (which is the driving force for the formation of FFLO state) d|^j»jji^es 



the physics. Indeed, some experimental indications of the existence of the FFLO state have been reportedav 
All of these materials have been argued to be unconventional superconductors and in this way differ from the case 
originally considered by Fulde, Ferrell, Larkin, and Ovchinnikov. Motivated by this possibility we have derived the 
Ginzburg-Landau (GL) free energy functional for unconventional superconductors in the presence of Zeeman splitting 
and impurity potential, and use it to study the possible FFLO phases. This is not a commonly used approach of 
examining the FFLO phase. To our knowledge it has previously only been discussed in the context of clean s-wave 
superconductors by Buzdin and KachkachiEZl However, as we show below, it represents a very powerful approach to 
study the FFLO phase since the simplicity of the resulting theory allows complexities such as non s-wave pairing, 
impurities, and even strong-coupling effects (though this is not done here) to be included. The stability of the various 
superconducting phases and to some degree, the topology of the superconducting phase diagram can be examined 
within this approach. These results help to clarify the nature of the FFLO phase and also can be used as a guide for 
a theory extended to all temperatures and magnetic fields. 

Using a functional integral formalism, we derive the GL free energy functional for single component singlet super- 
conductors in the presence of impurities and Zeeman fields in the weak coupling limit. The resultant GL free energy 
is valid near the second order normal to superconductor phase transition line in the (T, H) plane. This line will be 
denoted by [T, H(T)). The resulting instability to the FFLO phase appears readily within this approach due to the 
change of sign of the gradient term kIV^I 2 along the line [T, H(T)]. The point at which the coefficient k changes sign 
[denoted [T* , H(T*)]] is a tricritical point. At this point the normal, uniform superconducting, and FFLO phases all 
meet (see Fig. [l]). An intriguing feature of the weak-coupling clean limit is that the fourth order uniform term /3|'I , | 4 
also changes sign at the tricritical point. It is this term that determines the form of the FFLO phase. In particular, 
for a given momentum q both solutions ^ ~ e iq r and \& ~ e~ lqr are degenerate superconducting states at the normal 
to superconductor instability. The fourth order term breaks this degeneracy and selects either a cos(q • r) (LO phase) 
or a e lq r (FF phase) type order parameter. Since the magnitude of the order parameter is spatially uniform for the 
FF phase and vanishes at lines in real space for a LO phase, a negative (3 stabilizes the LO phase. Note that the 
complete GL free energy in this case requires inclusion of terms of the form |\&| 6 and |^| 2 |VW| 2 to be bounded. These 
considerations have appeared in the work of Buzdin and KachkachiU for conventional s-wave superconductors in the 
clean limit and are shown here to remain true for unconventional superconductors. 
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We further extend these considerations to include the effect of impurities. It is found that impurities suppress 
the FFLO phase for both conventional and unconventional superconductors. However, we also find that impurities 
lead to qualitatively different (TVff) phase diagrams for conventional (s-wave) and unconventional (non s-wave) 
superconductors (see Figs. ||and |J). This difference is most easily understood by looking at the coefficients n and (3. 
For conventional superconductors it is known that (3 is unchanged by non-magnetic impurities (this is a consequence 
of Anderson's theorem) while k is changed. The point on the [T, H(T)] line that n changes sign is pushed to lower 
temperatures with impurities (which illustrates that impurities suppress the FFLO phase) . Consequently, the normal 
to FFLO transition in the clean limit is replaced by a first order normal to uniform superconducting transition. For 
unconventional superconductors impurities change both the k and (3 coefficients, due to the absence of Anderson's 
theorem. It is found that the points on the [T, H (T)] line that k = and (3 = move to lower temperatures with 
increasing impurity concentration. The point [3 — is more rapidly suppressed than the point k = 0. This implies 
and that the initial instability into the FFLO phase is into a FF phase (\P ~ e* q r ) as opposed to the LO phase 
(if? ~ cos q • r) that is typically encountered. We are not aware of any other report of the stability of a FF phase. The 
results here also indicate that the first order normal to uniform superconductor phase transition does not occur for 
unconventional superconductors. 



II. GINZBURG-LANDAU THEORY 



Consider the Hamiltonian 

H = f dx^ft(x)[T(x) + C/(x) + 2<7//fl]tf (T (x) + V mt , (I) 

J a 

where T(x) represents the kinetic energy and takes the form — V 2 /(2to) — tF for a free electron and more generally 
takes the form T(x) = e(k = zV) for a band with dispersion e(k) measured from the Fermi energy tp , t/(x) is the 
disorder potential and satisfies (U(x)) — and (U(x)U(x')) = mWS d (x — x'), n; is the concentration of impurities, 
and /i = <7/j,b/2 is the magnetic moment of the electron. We will primarily be interested in singlet pairing in the 
ground state; so we neglect interactions between electrons with the same spin. The pairing interaction (Vi n t) is taken 
to have the separable form 

Vint = -V /k/k4+H, T C t _ k+fa C t _ k , + f d cl +§ >T (2) 

k,k',q 

f(k) describes the gap dependence on the Fermi surface and is defined to satisfy J2k l/k| 2 = 1- It is a l so understood 
that a cutoff exists in momentum space so that only electrons that are close enough to the Fermi surface interact with 
each other. After taking the appropriate Fourier transforms the interaction in real space takes the form 

V mt = -V J dxdx'dR/(x)r(x')*}(R+|)* t i (R-|)* i (R-y)* T (R+y) (3) 

where /(x) = -^X)k/ ke ' k X ' ^ ^ s ^ ne vonime of the system (note this definition of /(x) implies J <ix|/(x)| 2 = I). 
One may also describe the system using an Euclidean action in terms of Grassman variables: 

- Jj drVo J dxdx'dR/(x)/*(x')* T (R + f , r)* x (R - | ,r)* X (R - r)* T (R + £,r) (4) 
where Sq is the action for free electrons, and r is the imaginary time. The partition function is 

Z= f D^D^e-^'^. (5) 



We now decouple the quartic term in S by introducing a pair of Hubbard-Stratonovich fields A(R, r) and A(R, r), 
which will become the superconducting order parameter: 

5[*,*,A,A] = S - J^dr J rfRdr[A(R,T)/(r)* T (R+^,T)* i (R-^,r) 

+ A(R,r)/*(r)v|/ i (R-^r)vI/ T (R + |,r)]+^ ,3 dr J dR l -^^l. (6) 
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With this decoupling, the fermionic action becomes quadratic, and can be integrated out, after which we obtain an 
effective action in terms of the order parameter A(R, t): 



— r 13 lAfR t)\ 2 — 

S e {A,A} = / drdR 1 \; n -logZ[A,A], (7) 

where 

Z[A, A] = J L»¥D* 

e -S [*,^+J^dT / ( iRdr[A(R,r)/(r)* T (R+£,r)* i (R-£,r)+A(R,r)/*(r)* i (R-§,r)* T (R+§,r)] , g . 

The mean- field solution corresponds to the saddle point of S e [A, A]: 
SS e [A,A}. 

_ A 3 (R) 3 log Z [A, A] | 

= - J */(r)(^(R - §)* T (R + §)) As = 0, (9) 

where ()a, stands for quantum and thermal averaging in the presence of the pairing field A S (R). Here we have 
assumed a static saddle point so that A s has no r dependence. 

The functional integral formalism can be used to derive the effective Ginzburg-Landau free energy, in the vicinity of 
the second order normal to uniform superconductor transition line [T, H(T)]. This has been done for the short-range 
attractive interactions (that gives rise to s-wave pairing)cil. Our starting point is the effective action, Eq. (Q). Near 
[T,H(T)] ; we may make two simplifications: i) We may neglect the r dependence of A as we expect the thermal 
fluctuations to dominate the quantum fluctuations; ii) We may expand S e in powers of A. The quadratic terms take 
the form 



^ 2) [A, A] =/3 J dR |A ^ )|2 - J dRdR'Q(R, R')A(R) A(R') 



(10) 



where 



Q (R)R') = 6H ° S J |a-o 
^ V ' ' 0A(R)<SA(R') 



= J dnj dr 2 |drdr'/(r)r(r')(* i (R-^,T 2 )vI/ T (R + £,T 2 )vI/ T (R' + ^,r 1 )* i (R'-^,T 1 )} c 
= jf" dn J* dr 2 J drrfr7(r)/*(r')G 04 (R - ~ R' - ~; r 2 - n)Go, T (R + |, R + ^; r 2 - n) 
= J2 I drdr'f(r)t(r')G od (R-^,-R' - L ; ^„)G , T (R+ ^,R' + |; -iw n ). (11) 

The quartic term takes the form 

S< 4) [A,A] = ~ J dR 1 dR 2 dR 3 dR4i?(Ri,R 2 ,R3,R4)A(R 1 )A(R 2 )A(R 3 )A(R 4 ) (12) 

where 

fl(Ri,R 2 ,R 3 ,R4) = Jo dT 1 dr 2 dTzdT i J dridr 2 dr 3 dr4/(ri)/* (r 2 )/(r 3 )/* (r 4 ) 

G , T (R 4 - ^,Ri - ^-;r 4 - r 1 )G ,i(R 4 + T' R 3 + f J r 4 ~ T 3) 

G , T (R 2 - f ,R 3 - f ; r 2 - t 3 )Go,|(R2 + f , Ri + f ; r 2 - n). (13) 

We will need the sixth order term as well, whose explicit expression (that involves the product of six Green's functions) 
is not included here. The above quadratic and quartic terms apply for a particular impurity configuration. We 
will average over impurity distributions when calculating the form of the free energy. We assume that the gap 
function that appears above corresponds to the gap function averaged over impurities and that we can ignore impurity 
induced correlations in the gap function and between the gap function and the single particle Greens functions. To 
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proceed further the impurity averaged correlation functions (GG) and (GGGG) must be calculated. We determine 
these within theJiorn approximation and much of the derivation follows that of Werthamer for conventional s-wave 
superconductorsEi 



The impurity averaged normal Greens functions are 

Go,cr(k;iw„) = 



(14) 



where 



T = irmWN(0) and er is 1/2 (—1/2) for f {[). Consider the average Q(xi, yi; X2, y2? ioj n ) 



(Go,j,(xi,yi;iw„)G ! o,t(x2,y2; -iu> n ))imp rnO^i - y2,yi - y2,x 2 - y2\iu n ) due to translational invariance. Sum- 
ming the usual ladder diagrams shown inEa gives the self-consistent solution 



Q(Ri, R 2 , R3; iuj n ) — Go,x(Ri — R-2; «w n )Go ) f(R3; — «w„) 

+ rnW JdRG 0: |(Ri - R; iuj n )Q(0,R 2 - R,R 3 - R; zw n )G , T (R; -iu n ). 

This can be solved after taking the Fourier transforms with respect to Ri , R 2 , and R3 

G(ki,k2,k 3 ;iw„) = Go,t(ki; iu n )Go,i(k.2; -iu n ) 

^ki,-k 2 



(15) 



rnWGo,t (ks:iuj n )Go,i (ki+k2+k3; — iw„ ) 
l- Z T? 1 Xi k ^ 'T(k;i^)Goa(k+k 2 +k3;-iw„) 



Substituting this result into Eq. 11 gives 
Q(R,R) = 



V Ec 



iq(R-R' 



E 



G , T (k- 



l )G ,|(k'-§;^„) 



VS ] 



k,k' 



k,k',M„ /k/k' 

n i WGo, 1 (k'+%:iiO n )G ,i(k-%;-tw n ) 

■Hr- E P g o,t(p-| ; iw ") G o,i(p+f 



(16) 



(17) 



Note that the form of the vertex corrections found here is not the same as for conventional s-wave superconduc- 
tors. In particular, when deriving the terms up to second order in the gradients, the vertex corrections vanish for 
unconventional superconductors. However, for higher order gradient terms, the vertex corrections are not zero. 

The impurity averaged correlation function that appears in SS4 is less straightforward to calculate. For the terms in 
the free energy that are fourth order in the order parameter we consider the only up to second order in the gradients 
of thc,order parameter. In this case the non-zero diagrams have the same form are those that contribute in the s-wave 
caseEl 

After performing the appropriate Taylor series expansions the following GL free energy for hexagonal and square 
lattices is found [this expression is valid only for non s-wave superconductors ((/k) = 0)] 



F = a\A\ 2 + P\A\ 4 + k|VA| 2 + £|V 2 A| 2 + M |A| 2 |VA| 2 + ?7[(A*) 2 (VA) 2 + (A) 2 (VA*) 2 ] 



the coefficients arc 



HA| 6 + ^|(V 2 -V 2 )A| 2 



a = -N{0) [ln(T c °/T) + nKi - txK x (r = 0, B = 0)] , 



(18) 
(19) 



/?= ^V(0) ({|/(k) |4 ) ^ 3 _ r ^ 4)) (20) 



^V(0)(vi(k)|/(k)| 2 ) 
K = i 3 K 3, (21) 



7TiV(0)(|/(k)| 2 vi(k)) 

5 = ^ K 5 , (22) 

7rjv(0)(vj(k)|/(k)| 4 ) ^ r 
/x = 8r?= (A 5 - (|/(k)|4) Jf 6 ), (23) 
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6 = v - (K 5 + TK 6 ), 



(24) 



K 6 + 2T 2 K 7 ), 



(25) 



where 



K n = (2T) 1 - n —Re ( V . 1 . 



(26) 



lattice. For an orthorhombic lattice the following terms also appear in the free energy 

SF = *(|V X A| 2 - |V y A| 2 ) + /2|A| 2 (|V^A| 2 - |V,A| 2 ) 

+ r?[((V x A) 2 - (V,^) 2 )(A*) 2 + ((V.A*) 2 - (V,A*) 2 )(A) 2 ] 
+ h {[(V 2 + V 2 )A][(V 2 - V 2 )A]* + [(V 2 + V 2 )A]*[(V 2 - V 2 )A]} 



The 8 term does not appear for a hexagonal 



(27) 



The coefficients k, jj, and i) are given by k, /i, ij, with v 2 ^ replaced by (u 2 — v 2 ) respectively and the coefficient 62 is 



given by S with Vj_ replaced by vj_ (v x 



The free energy is the main result of this paper. 



III. 5"- WAVE SUPERCONDUCTORS 



Here we review the known results about the FFLO phase for a conventional s-wave superconductor and show how it 
arises from the free energy. As mentioned in the Introduction the only two coefficients that are required to study the 
instability from the normal phase into the FFLO phase near the tricritical point are k and (3. For a conventional s-waiie 
superconductor (for which /k is a constant) these can easily be determined by following Werthamer's derivation! 2 ^ 
(note that the above free energy does not apply here since (/k) 7^ 0): 



7rJV(0)(vj(k)) - 

n K 3 



(28) 



ttN(O) , 



(29) 



where 



K 3 = {2T)- 2 ^Re 



OO 



^ 2 l 2irT> ^ ' 2 1 2ttT ~ 2ttT 



(30) 



The coefficient (3 does not depend upon the impurity concentration in agreement with Anderson's theorem. The 
second order normal to uniform superconductor phase line [T, H(T)] is given by a(r = 0) = and is shown in Fig. |J. 
Note that once k < or < the phase line [T,H(T)] no longer denotes the true normal to superconductor phase 
line. Numerical evaluation of k s and (3 show that in the clean limit both k and [3 vanish at T = 0.56T C . This point is 
the tricritical point. The phase diagram is shown in Fig. El When impurities are added numerical evaluations show 
k vanishes at a lower temperature than (3 vanishes (note that the presence of impurities does not change [3). In this 
case there is a first order superconducting transition to a homogeneous phase for T < 0.56T C . Once (3 < then the 
normal to superconductor instability line is no longer given by Fig. |[ To determine whether there exists an FFLO 
phase for some arbitrary impurity concentration requires a calculation that goes beyond the GL free energy presented 
here. This is because the G.L. theory is only valid close to the transition at T = 0.56T C (note that the G.L. theory can 
be used to study the transition into the FFLO phase if the impurity concentration is small enough) . The calculations 
of Bulaevskii and Guseinov for layered superconductors indicate that there exists no FFLO phase at T = when 
r/T c > 0.6l3. The qualitative phase diagram in the presence of impurities is shown in Fig. |l 
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IV. UNCONVENTIONAL SUPERCONDUCTORS 



The last Section demonstrated that for T near 0.56T C the G.L. theory accurately reproduced the phase diagram for 
conventional superconductors. Here we apply the same approach to unconventional superconductors where it turns 
out the G.L. theory is more powerful. This arises because the transition from the normal state to the superconducting 
state is second order for all fields and impurity concentrations for unconventional superconductors. For conventional 
superconductors this transition is sometimes first order, which limits the applicability of G.L. theory. For example, 
the G.L. theory for unconventional superconductors can give the maximum impurity concentration that allows the 
FFLO phase to exist; it was not able to do this for s-wave superconductors. As a concrete example we study a d-wave 
superconductor (/k °c k^.~ky) with a cylindrical Fermi surface. Choosing some other /k will not change the qualitative 
form of the phase diagrams (provided (/k) = 0). The clean limit phase diagram is qualitatively the same as that for 
the s-wave case. In fact, the clean limit theory indicates that the FFLO phase appears for T < 0.56T C independent 
of the order parameter symmetry. When impurities are added the main conclusion is that k vanishes at a higher 
temperature than (3 does when these quantities are evaluated on the phase boundary [T, H(T)] (see Fig. || and Fig. ^ 
for the temperature evolution of k and (3) . This implies that there is no first order transition from the normal state to 
a uniform superconducting state, but rather a second order transition into a FFLO phase. The temperature at which 
k = gives the maximum temperature that allows the existence of the FFLO phase (we call this temperature Tp). 
This is in sharp contrast to what happens in s-wave superconductors, where a first-order phase boundary separates 
the normal and uniform superconducting state for temperatures just above the point at which (3 vanishes. Fig. |?] 
shows how Tp varies as the impurity concentration is increased (note that impurities also suppress T c , hence we plot 
Tp/T c as a function of T c /T c q since the latter is an experimentally measurable quantity). This figure indicates that 
the FFLO phase survives a considerable impurity concentration, only for T/T c > 0.6 does the FFLO phase cease to 
exist (superconductivity is destroyed when T/T c > 0.88). 

The structure of the FFLO phase can also be addressed within the G.L. theory in the neighborhood of tricritical 
point (given by k — on the phase line [T, H(T)]). To do this we compare the free energy for three different phases: 
(1) Ai = A (uniform phase), (2) A 2 = A e iqr (FF phase), and (3)A 3 = A cos(q-r) (LO phase). We note here that 
terms 0(A 6 ) do not need to be included in the calculations done here because the terms 0(A 4 ) are positive. Also 
the minimization with respect to the orientation of q for A 2 and A3 implies that q is oriented along the nodes fpc all 
impurity concentrations. This agrees with earlier calculations that go beyond the GL theory in the clean limitli3'[L3. 
After minimizing with respect to q the resulting free energies are 

Fi= a|A| 2 +/?|A| 4 (31) 

^Wa-^)|A| 2 +/? 2 |A| 4 (32) 
^ = (a-^)|A| 2 +/3 3 |A| 4 (33) 

(34) 

where (3 2 = (3 - and (3 3 = |/3 - ^f 2 (where we have used (i = 877). It is clear that when n = 0, Fi = F 2 < F 3 
which implies that a — and k = gives the tricritical point where the normal, uniform superconducting, and FF 
superconducting phases meet. Note that the FF phase is stable while the LO phase is not at the tricritical point 

since (3 > (the FF phase is never stable in the s-wave case). For temperatures below that of the tricritical point 

2 

the phase transition from the normal phase into the FFLO phase is given by a — jr. Intriguingly, along this phase 
line (3 changes sign. This implies a first order transition between the FF and the LO phase. This can be seen by 
comparing F2 and Fa which is equivalent to comparing /3 2 and (3 3 . If [3 = and k < then < (3 3 < /? 2 which implies 
F3 < F2 implying that there is a phase transition between the FF and the LO phases. For the singlet superconductors 
considered here the FF phase should exhibit a spin current. 

A detailed calculation was carried out for an impurity concentration for which T c /T c q = 0.573 (where T c q is the 
transition temperature with no impurities present). Fig. || shows the resulting phase diagram calculated within the 
G.L. theory. Note that G.L. theory in this case gives a reasonable description of the phase diagram in this case because 
it is valid along the entire normal to superconducting phase transition boundary. For other impurity concentrations 
the phase diagram is similar. If multiple scattering from the impurities becomes important (e.g. going beyond the 
Born approximation to a T=matrix treatment), then it is found that the region of the phase diagram where the FF 
phase appears is decreasednll. 



G 



V. CONCLUSIONS 



In conclusion we have derived the GL free energy for singlet superconductors in the presence of a Zeeman field 
and non-magnetic impurities. This free energy was used to examine the resulting phase diagram. It was shown that 
the phase diagrams for unconventional superconductors and conventional superconductors are qualitatively different 
in the presence of impurities. In particular the first order normal to uniform superconductor phase transition that 
exists for conventional superconductors does not exist for unconventional superconductors. Also, for unconventional 
superconductors impurities induce a change in the structure of the FFLO phase. In the clean limit the FFLO phase is 
described by an order parameter of the form cos(qr) (LO) while impurities stabilize a e lqr (FF) type order parameter. 
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Magnetic 
Field 




Temperature 

FIG. 1. Qualitative phase diagram for clean superconductors. The cosqr form of the order parameter is only valid near the 
normal to FFLO transition line. The direction of q in the FFLO phase may also depend upon temperature (see for example 
Ref. [18]). 
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Temperature 

FIG. 2. Qualitative phase diagram for conventional superconductors with non-magnetic impurities. 
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FIG. 3. Qualitative phase diagram for unconventional superconductors with non-magnetic impurities. 
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' 0.0 0.2 0.4 0.6 0.8 1.0 

FIG. 4. Second order normal to uniform superconducting phase boundary for s-wave superconductors. This line defines 
[T, H (T)] when there are no impurities present. This is the boundary on which k and /3 are determined. Note that for 
T < 0.56T C this phase boundary will not coincide with the actual normal to superconducting phase boundary (the transition 
will either be to a non-uniform (FFLO) phase or will be first order). 




FIG. 5. k(T)/k(T c ) on the phase line [T,H(T)] for d-wave superconductors. The curves from top to bottom at T = T c 
correspond to T/T c = 0.6, T/T c = 0.3, and T/T c = 0.0 respectively. 
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FIG. 6. /3(T)/(3(T C ) on the phase line [T,H(T)] for d-wave superconductors. The curves from top to bottom at T — T c 
correspond to T/T c = 0.6, T/T c = 0.3, and T/T c = 0.0 respectively. 



FIG. 7. The maximum temperature (Tp) for which the FFLO phase can exist as a function of the transition temperature 
(which is suppressed by impurities). 




FIG. 8. Phase diagram for a d-wave superconductor with T c suppressed by impurities such that T c /T c o = 0.573. FF refers 
to an order parameter of the form A = Aoe lqr while LO refers to an order parameter of the form A = Ao cos(q • r). 
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